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The history and recent progress of atmospheric neutrinos are reviewed. An emphasis is placed on results from 
experiments other than Super-Kamiokande. 
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The Super-Kamiokande experiment [0 has used 
their measurements of atmospheric neutrinos to 
persuasively and clearly demonstrate the ex- 
istence of neutrino oscillations, with excellent 
agreement with the hypothesis i^^ Vr- With 
lower statistics, other experiments have measured 
atmospheric neutrinos for a long time and con- 
tinue to do so. In this paper, after reviewing the 
creation and measurement of atmospheric neutri- 
nos, I discuss the history of atmospheric neutrinos 
along with two possible alternative scenarios. I 
then review the most recent results from Baksan, 
Soudan 2 and MACRO. 

2. The Creation of Atmospheric Neutrinos 

Atmospheric neutrinos originate from the de- 
cays of tt's, K's and ^'s produced when cosmic 
rays hit the atmosphere and interact. There are 
also a smaller number of neutrinos produced by 
the decay of charmed particles, r's and other high 
mass particles, but their detection has not been 
demonstrated. K and tt decay give mainly z/^'s 
while iJ, decay gives both i/^'s and ly'^s. The ^'s 
themselves come from tt and K decay, so at low 
energy (below 2 GeV where all of the muons de- 
cay before they hit the earth), the flux of each 
flavor neutrino (jfi^iJ. + ^ij.) ■ 4>iye + f'e) occurs in 
the ratio 2:1. The experimental value for this 
ratio (in the absence of oscillations) for charged 
current interactions will vary from 2 for several 
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1. Containment differences between and e 

2. charged current threshold 

3. and differences due to their different 
cross sections and the fact that ratio 
is about 1.2 

4. Some additional v^s from K'^^ and K'^^ de- 
cays 

5. Higher energy /i's hit the earth and lose 
most of their energy before decaying. 

For detectors built to study proton decay, the 
first factor is the most important, and the ratio 
is under 2. For high energy neutrino telescopes, 
the last factor causes the ratio to increase to large 
values. 

The density profile of the atmosphere affects 
the geometry of the source. The overburden, and 
hence the pressure, [|| goes as 

X =Xq cxp i-pgh/kT) = Xa exp {-h/ho) (1) 

where the latter expression is used to define the 
scale height ho, which is about 8.4 km at sea 
level. Since T and hence ho depend on altitude, 
ho decreases to 6.4 km near the tropopause, where 
many atmospheric v are created. Also, the air 
temperature induces seasonal variations of muon 
and neutrino fluxes. Muon seasonal variations 
are typically ±2%, but since high energy neutri- 
nos come primarily from if's, this effect is small 
for high energy neutrinos. A bigger effect at low 
energy is the solar cycle variation, which occurs 



because the magnetic field of the solar wind pre- 
vents many low energy cosmic rays from reaching 
the earth's location in the solar system. 

If all neutrinos came from the same height in 
the atmosphere, there would be a strict relation- 
ship between the local neutrino zenith angle in 
an underground detector and the neutrino path- 
length L which is relevant for neutrino oscillation 
analysis. Those production height distributions 
vary with both neutrino angle and energy. A 
useful parameterization has been performed by 
Ruddickl^. A comparison of his parameteriza- 
tion with a full Monte Carlo at two choices of 
energy and angle are shown in Figure 1. 



3. Detection of Atmospheric v 

Two classes of atmospheric neutrino events are 
considered in this paper. A neutrino may inter- 
act in an underground detector, or it may interact 
outside the detector and make a muon which is 
detected. Events of the first category include the 
contained events and partially contained events 
of Soudan 2, 1MB and Kamiokande and the semi- 
contained events of MACRO. Events of the sec- 
ond category include the upward throughgoing 
/i's and stopping /i's of MACRO and Kamiokande 
and the horizontal /i's of KGF, Soudan 2 and Fre- 
jus. An attempt to catalog the world total of at- 
mospheric z^'s is given in Table 1. Only neutrino 
candidates which are background free or back- 
ground subtracted have been included. It is clear 
that Super-Kamiokande dominates the total. 
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Table 1 

Atmospheric neutrino event totals. The *'s rep- 
resent experiments that are still taking data. The 



total, not including Super-Kamiokande is 6202 v 
events 



Figure 1. Production height distributions for two 

neutrino energies and angles. 4. History 



The first reported atmospheric v was measured 
in the KGF experiment H using a set of telescope 



counters consisting of iron, flash tubes and scin- 
tillator. That telescope was operated on several 
levels of the KGF mine, but started at 7600 MWE 
level where the nucleon decay and monopole de- 
tectors were later built. The measured angu- 
lar distribution in the latter detectors shows a 
clear separation between atmospheric /x's and v 
induced /i's. This is shown in Figure ||. Even- 
tually, 100 contained events and 229 v induced 
horizontal /i's were measured at KGF. 
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muon and neutrino induced muon can be clearly 
separated was about 50°. 
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Figure 3. Koshiba's slide from 1988 Fermilab 
Workshop. 



Figure 2. Zenith angle distribution of muons in 
KGF. 

The first recorded atmospheric v was measured 
on 23 February, 1965 in the South African Neu- 
trino Detector (CWI/SAND) built by the Case- 
Western- Irvine Group. Q Again, the signature for 
an atmospheric v was the projected zenith angle. 
Due to the extreme depth of that detector (8890 
MWE), the zenith angle where an atmospheric 



More atmospheric neutrino data became avail- 
able in the 1980's from large detectors which were 
built to search for nucleon decay. The first hint 
of the atmospheric neutrino problem, which was 
also known as the "too few nu mu" problem came 
in the IMBl data. They measured the number 
of delayed coincidences due to muon decays. H 
This distinguished v^Jl — s- /i^p and P^p — > iTn 
followed by /X — > evv from v^n —>■ e~p and 
v^p — >• e+n events. They reported, "The simu- 
lation predicts 34% ± 1% of the events should 
have an identified muon decay while our data has 
26% ± 3%. This discrepancy could be a statis- 
tical fluctuation or a systematic error due to (i) 
an incorrect assumption as to the ratio of muon 
i/'s to electron z/'s in the atmospheric fluxes, (ii) 
an incorrect estimate of the efficiency for our ob- 
serving a muon decay, or (iii) some other as-yet- 



unaccounted-for physics." It has turned out to be 
the latter. 

To learn from an anomaly, there needs to be not 
only results inconsistent with expectation, but 
also someone to take it seriously and to carefully 
study alternative explanations. The first detailed 
evaluation of the atmospheric neutrino problem 
that I personally saw was a discussion by Koshiba 
at the 1988 New Directions in Neutrino Physics 
Fermilab Workshop. There he showed the trans- 
parency in Figure ^ He compared the "ratio-of- 
ratios" defined as 



for several experiments and several different 
methods. The Water Cerenkov Detectors could 
measure the ratio two different ways, with muon 
decays and with ring fits. The initial measure- 
ments of the iron calorimeters, Frejus and NU- 
SEX got higher values consistent with unity, but 
NUSEX had quite low statistics, and even Frejus 
had a noticeable i/^ deficit if one looked only at 
the fully contained events. The period 1988-1995 
was the period of the "ratio-of-ratios" or the "at- 
mospheric neutrino anomaly" . Several attempts 
to understand the anomaly either as a systematic 
effect or as an error in the atmospheric neutrino 
Monte Carlos did not succeed. 

In 1994, the Kamiokande experiment pO| 
showed another data set, their Multi-GeV data. 
In the Multi-GeV data, not only was there a 
deficit of v^, and a value of R below unity, but 
there was also a zenith angle distribution consis- 
tent with neutrino oscillations. By contrast, their 
zenith angular distribution for the Sub-GeV con- 
tained events was flat. This suggested a higher 
value of Ato^ than currently seen. It has been 
widely noticed that the parameter space plots ( 
in sin^(26') and Am^) for Kamiokande and Super- 
Kamiokande did not overlap. It is interesting to 
compare the zenith angle plots of Kamiokande 
data and Super- Kamiokande data for Sub-GeV 
and Multi-GeV, e and ji. In such a compari- 
son, no data points disagree by much more than 
one sigma. The disagreement in parameter space 
comes about in the fits. It is worthwhile to point 
out that fits for Am^ are not gaussian, and that 



it is possible to have multiple solutions similar to 
the degeneracies now facing experiments planning 
to measure ^13 and 5c p- 

5. Gedanken History 

Let's imagine the status of atmospheric neutri- 
nos for two scenarios in which history had been 
different. 

The largest underground experiments for nu- 
cleon decay and magnetic monopole detection 
were motivated by Grand Unified Theories. In 
the absence of such motivation, the only exper- 
iments which measured atmospheric v^s would 
have been CWI/SAND, Baksan and LVD. So- 
lar neutrino experiments would have proceeded 
on their same time scale (or even faster). There 
would likely have been a greater interest in the 
reported neutrino oscillation signature from the 
LSND experiment [r^, as well as in proposed 
short and intermediate baseline experiments such 
as BNL to Long Island, CERN to JURA and the 
never-realized Fermilab COSMOS experiment. 
The solar neutrino results in the 80's and 90's 
might have motivated some forward thinking in- 
dividuals to propose a large underground exper- 
iment to measure atmospheric neutrinos. There 
would probably have not been great enthusiasm 
for such an experiment, most sensitive to large 
mixing angles, until after the SNO results in 2002. 
Even then, I think it would have been a hard 
project to realize. 

History would also have been quite different 
if the Super-Kamiokande accident had happened 
during its first fill in 1995. Confidence in rebuild- 
ing the detector might not have been possible in 
the absence of the 5 successful years of running 
the detector. The K2K run would not have taken 
place. MINOS, which had already been approved, 
would have continued to plan with its high energy 
beam, most sensitive to larger values of Am^. 
New data from MACRO and Soudan 2 would be 
tending to support lower values of Am^ , but given 
the non-gaussian (and hence unintuitive) nature 
of Am^ fits, the situation would be fairly con- 
fused, which would have lead skeptics to doubt 
conclusions about atmospheric neutrino oscilla- 
tions. The latest Soudan 2 and MACRO data 



analyses would be quite relevant in trying to sort 
out the situation. 



6. Recent Results from Baksan 



The Baksan Underground Scintillator Tele- 
scope has been taking data since the 1980's, with 
its four layers of scintillator detectors under a 
mountain at a minimum 850MWE. At Neutrino 
2000, they reported on data from December 1978 
to January 2000, corresponding to 15.7 years of 
livetime. They measured 801 upward /^'s with an 
expected rate (in the absence of oscillations) of 
941.6. The zenith angle distribution was not in 
close agreement with oscillation fits, though the 
rates agree with other experiments. They now 
have about another 10% increase in statistics, but 
will not present a new analysis until next year or 
later, when the increase is 25% or more. 

7. Recent Results from Soudan 2 

Soudan 2 is a very fine-grained drifting 
calorimeter with drift cells located in a honey- 
comb pattern of iron plates. Since Neutrino 2000, 
the Soudan 2 experiment has finished taking data. 
It finished in July 2001 with a total 5.91 fidu- 
cial kt-year for contained event data. Details of 
the current status of Soudan 2 analyses for atmo- 
spheric neutrinos can be found in Reference [|l8j . 
A recent feature of the Soudan 2 data analysis is 
that a 15% problem with the electron energy scale 
has been resolved. Data analysis now includes 
the partially contained events, and parameter 
space fitting is done with a Feldman-Cousins [p^ 
type of analysis. A feature of Soudan 2 which 
has been exploited for some time is the ability 
to measure recoil protons, which appear in this 
very fine-grained detector as short, straight, heav- 
ily ionizing tracks from the main vertex. The 
events which have an identified recoil proton al- 
low the neutrino zenith angle, and hence the L 
relevant for L/E analyses, to be reconstructed 
with much greater accuracy. Soudan 2 defines 
its high-resolution sample to be the high energy 
quasi-elastics, the low energy quasielastics with a 
recoil proton and the high energy multiprongs. 
The zenith angle distributions of the electrons 
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Figure 4. Zenith angle distributions in Soudan 2 
for c (top) and /i (bottom). Shown are data 
points, and histograms with and without oscil- 
lation fits. A small background is also shown. 



and muons are shown in Figure I and the L/E 
distribution in Figure ||. 




0.5 1 1,5 2 2.5 3 5.5 4 4.5 

L/E, rM.jor,s. 



Figure 5. log^o L/E distribution in Soudan 2 for 
e(top) and /i(bottom). The histograms have the 
same meaning as in Figure 4. 

The Feldman- Cousins analysis involves calcu- 
lating the likelihood of the data as function of 
sin^(26') and Am^ parameters for the data and 
for Monte Carlo data. A confidence level diagram 
is calculated from the difference in log-likelihood. 
This is shown in Figure H The best fit occurs at 
a value of Am^ = o.OlO eV'^ and sit? {29) = 0.97. 
ft should be noted that the valley near the best 
fit is quite flat, particularly towards lower values 
of Ato^. 

8. Recent Results from MACRO 

The MACRO experiment is a large area detec- 
tor in the Gran Sasso laboratory, consisting of 3 
towers of scintillation counter and 14 horizontal 




Figure 6. Likelihood surfaces for Soudan 2 data. 



planes of streamer tubes. MACRO completed its 
data taking in the last two years, finishing with 
5.52 live years in December 2000. With excel- 
lent timing, MACRO was able to separate up- 
ward and downward going events, but the trig- 
ger gave poor acceptance for horizontal muons. 
The categories of events in MACRO were 809 -|- 
54 [background subtracted H- background] (1122 
[Monte Carlo]) upthroughgoing muons, 154 -1-7 
(285) internal upgoing (lU) events, and 262-MO 
(375) upgoing stopping muons (UGS) plus inter- 
nal downgoing muons (ID). 

The shapes of the zenith angle distributions of 
the MACRO event samples are sensitive to neu- 
trino oscillations. For the throughgoing muons, 
the ratio of events with —1 < cos6' < —0.7 to 
events with —0.4 < cos6' < can be used to dis- 
tinguish Vp, — > Vr from jAj Vs oscillations be- 
cause of a matter effect . The has a neutral 
current interaction in matter while the Vs does 
not, and this difference in interaction leads to a 



different expected oscillation probability. Based 
on the ratio test, Vs oscillation has only 

a 0.033% probability of fitting the data, and is 
disfavored by > 99%CL compared to the best fit 
i/^ — > v-r oscillation. 
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resolution using the limited streamer tubes in 
drift mode, and through an analysis of the multi- 
ple scattering, obtained an energy estimate on a 
muon by muon basis. With that energy, and the 
zenith angle which gives an estimate of the neu- 
trino flight path, they have made an L jE^, distri- 
bution. This is shown in Figure ^ along with the 
expected distribution for the best fit oscillation 
parameters. 

The allowed parameters (Am^, sin^(20)) ac- 
cording to the low energy and high energy 
MACRO analyses are shown in Figure || along 
with the analysis of Soudan 2. The results of 
both experiments agree with each other and with 
Super-Kamiokande. Each experiment also rejects 
the null hypothesis of no oscillations with high 
probability. 



Figure 7. Angular Distributions of MACRO low 
energy neutrino data, along with expected shapes 
with and without oscillations. The shaded area 
shows the uncertainty of the flux estimate. 

The angular distribution of the MACRO low 
energy data is given in Figure ^. For both the 
lU and UGS-I-ID data sets, the data falls below 
the expectation (without oscillations), even tak- 
ing into account the normalization uncertainty of 
the expectation. And together with the through- 
going data, they form an acceptable fit to the 
same neutrino oscillation parameters. 

MACRO has also recently shown new results 
which use muon energy in the analysis. The 
muon energy is not measured calorimetrically or 
magnetically. However, the multiple scattering in 
the detector is dependent on the energy on aver- 
age in a known way. The projected displacement 
from a straight line fit for a relativistic muon is 
inversely proportional to the muon momentum. 
They have obtained an improvement of the space 



o 1 - 




oa - 



log„,L(km)/E„(GeV) 



Figure 8. LjE distribution of MACRO data 
based on multiple scattering estimation. 




Figure 9. Comparison of parameter space esti- 
mates for MACRO and Soudan 2 



9. Summary 

There are several experiments besides Super- 
Kamiokande which collectively have observed at- 
mospheric neutrinos. Analyses of these data 
agree with the essential Super-Kamiokande re- 
sult: neutrino oscillations are required to account 
for the data. Taken as a whole, contained event 
data from 1MB, Kamiokande, Soudan 2, and Fre- 
jus all agree with a 30%-40% deficit of i^^ in- 
duced events, and also agree with the higher 
statistics Super-Kamiokande results. New fea- 
tures from Soudan 2 data strengthens that con- 
clusion. Soudan 2 now has an up/down differ- 
ence in its contained events which is statistically 
significant, and supports the region of parame- 
ter space found by Super-Kamiokande. With its 
recoil proton identification and good angular res- 
olution, Soudan 2 has the resolution, but not the 
statistics, to see the "reappearance" in the L/E 
distribution. 

Detailed analyses of v induced muons support 
a similar conclusion. All experiments see a deficit 



Figure 10. Comparison of zenith angle distribu- 
tions of throughgoing muons 

of Vp, consistent with the Super-Kamiokande ob- 
servations, and an angular distribution which 
is more consistent with the neutrino oscilla- 
tion hypothesis than with the null hypothesis, 
as shown in Figure |l^, a comparison of data 
from Baksan, Kamiokande, MACRO and Super- 
Kamiokande. |^ It is interesting to point out that 
the probability of the best fit for all four exper- 
iments is less than 30% each, so it is worthwhile 
to continue to investigate whether some aspect of 
the physics is not being modeled correctly. 

It is left to another speaker at Neutrino 2002 
to consider the proposals for future atmospheric 
neutrino oscillation projects.]^ However, I want 
to make a general observation, that much of 
the progress in atmospheric neutrinos has come 



from experiments that were designed to search 
for GUT predicted nucleon decay. The motiva- 
tion for much larger experiments to continue that 
search is quite strong, and we as a field have the 
technical means to accomplish such a search in 
our lifetime. I hope that search takes place. If 
it does, a much firmer understanding of atmo- 
spheric neutrinos will be an inevitable outgrowth 
of such an effort. 
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